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A Unified Model of the Prompt Optical Emission of Gamma-Ray 

Bursts 

Hirotsugu Doi 1 , Kentaro Takami 1 , and Ryo Yamazaki 1 
ABSTRACT 

The observational diversity of optical emission, which coincides with prompt 
gamma-ray bursts (GRBs), has been discovered in the recent Swift era. We show 
that on the assumption of the synchrotron radiation for the observed energy 
range below the X-ray band, the observed diversity can be explained using the 
internal shock model by taking into account a high-latitude emission and the 
spectral change due to the synchrotron self-absorption. It may even be possible 
in our model to include bright optical flashes found, e.g., in GRB 990123. The 
prediction of our model is that the spectral index in the optical band is dependent 
on whether the optical light curve correlates with those in the X-rays and/or 7- 
rays or not, which will be tested in the near future. 

Subject headings: gamma rays: bursts — gamma rays: theory 
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Introduction 

Thanks to Swift ob servations, 7- ray bursts (GRBs) can be rapidly followed-up in various 



observation bands (see Izhanjl2007l for a recent 



review 



In particular, some e yents have 
been observed in the optical band during t he prompt 7-ray-active phase (e.g., lYost et al. 



20061 : iRvkoff et alibood : iRoming et allhood ) . There is observational diversity in the prompt 
optical emission and several categories can be found as described in the following: 

(i) The optical light curve synchronizes with the X-ray/7-ray light curve with the 



optic al-to-X / 7-rav flux ratio almost constant with time (e.g., GRB 041219a: IVestrand et al. 
2005h . The peak flux ratio R pcak = F^ ak /F^ k is from a few to ten. 
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(ii) The optical light curve is smooth and mu ch less variable compared with X-ray/7-ray 
light curve (e.g., GRB 060210; IStanek et al.ll2007l ). The peak flux ratio -R pea k is almost equal 
to or smaller than unity for GRB 060210. 

(iii) The optical light curve is a superposition of a smooth com ponent and a variable 
one synchronizing with X-ray/7-ray light curves (e.g., GRB 050820A; IVestrand et al.ll2006l ) 
. This may be an intermediate case between cases (i) and (ii). 

(iv) A few events sh ow bright optical fla shes just after the 7-ray-active pha se. A typical 
exam ple is GRB 990123 jAkerlof et al.lll999h . an d for this event R ? P ^ ~ 2 x 10 2 fariggs et al. 
1999h . GRB 060 11 IB may also be in this class fclotz et al.lbood ). 



(v) For a fraction of the events, UVOT on Swift detected no optical counterpart within 
IP 3 sec after the BA T trigger, and it gave strict upper limits on the prompt optical emission 
(IRoming et al.ll2006l ). Although it is difficult to discuss quantitatively, some events may have 
dim prompt optical emission. 



The origin of the prompt optical emission has been widely discussed (e.g. JPanaitescu fc Meszaros 



20071 : IPanaitescu fc Kumail 12006^ The popular inter pretation is 



20001 : iFan et al.ll2005l : [Wei 

that of case (i), which postulates an internal shock origin (IVestrand et ah] 12005), while case 



m 



1999; IWang et al.ll2000l ; IZhang et al.ll2003l : iNakar & Piranll2004 



iv) postulates an external reverse shock emission origin fe.g..lSari fe PiraJl 99*91: 



McMahon et al 



VIeszaros fc Rees 



20061 ). Case 



(iii) is a superposition of both components. In this Letter, however, we will show that all of 
the cases are explained in the internal shock model. 



2. Prompt emission model and optical/X-ray light curves 



Basically, we adopt the same model as in the previous works (jYamazaki et al.l 12004 



Toma et al.ll2005al lbl). However, in this Letter, a much simplified version is considered. The 



central engine launches A^-emitting shells in the same direction, with the Lorentz factor 
7 = (1 — /3 2 )~ 1//2 and the opening half-angle of A6* t ot- We introduce the spherical coordinate 
system (t, r, ip) in the central engine frame, where the origin is the location of the central 
engine and $ = is the axis of the whole jet. We assume that an observer sees the jet on-beam 
(i.e., in notation of our previous works, A6^ b = A6* tot , i? w = 0, and $ b s = 8^ = 0, then 
A0 (j) = 7r). The departure time of each emitting shell t^ p is assumed to be homogeneously 
random between t = and t = tdm, where tdm- is the duration time of the central engine in 
its own frame. The jet emission instantaneously arises at r = r$. Then the observed flux 
F U (T) at an observer time T and at an observed frequency v is a superposition of each shell 



- 3- 



emission calculated as (IToma et al.ll2005bl ) 

F^(T) = K r U) - 2 f(r U) u/ 1 ) 



(i; 



where 

r w = (ci 2 /r )(T-t£ p ) (2) 

and Kq is the normalization constant. Here T = is chosen as the time of arrival at the 
observer of a photon emitted at the origin at t = 0. In our model, synchrotron radiation is 
considered. Assuming the standard internal shock model and typical model parameters, the 
break frequencies at the on-beam observer are calculated as v m ~ 30 keV and v a ~ 10 eV, 
where we assume the bulk Lorentz factor V = 300, the Lorentz factor of the internal shock 
T s h = 5, the outflow luminosity L m = 1 x 10 52 erg s _1 , the observed typical variability 
timescale St — 1 sec, and the shock microphysics parameters e e = Eb = 0.5 (see eqs.(l)-(3) 



in 



Weill2007l ). The cooling frequency at the on-beam observer u c is found to be much smaller 



than v a . Note that considering several uncertainties of the parameters arising in the internal 
shock model, v a could be around the R-band frequency = 1.7 eV. Then the sp ectral shape 
function, f(u'), is given by (e.g., IPanaitescu fc Meszarosll2000l ; iFan fc Weill2004l ) 



WKf 12 

WJv'a)- 1,2 {v'/0- pl2 



y a < v' < v' m ) 



(3) 



which is different from that for the case v' a < v' c < v' m (ISari et al.l 119981 ). One should remark 
that for v 1 < v' a , has a steep slope because of the synchrotron self-absorption. We don't 
consider the spectrum in the range v' < v' c in Eq. (jHJ), because we assume v' c <C v' a based on 
the above order-of-magnitude estimation and then our result does not depend on the value 
of v' c . In order to more accurately reproduce the observational properties such as -R pe ak, we 
may have to introduce another emission component that is responsible for the high-energy 
X-rays a nd /or 7-rays in th e observer frame, whose spectrum is characterized by the Band 
function (IBand et al.lll993l ). However such a correction is neglected here for simplicity. The 
pulse starting and ending times of each shell emission are 



start 



rnU) 
end 



*£, + (ro//fc)(l-/J) , 
t ( £ p + (r /Pc)(l-(3cosA6 t , 



(4) 
(5) 



In this Letter, the cosmological effect is neglected (i.e., we set z = 0). In the following, we 
adopt the fiducial parameters r = 1 x 10 15 cm, 7 = 100, A9 to t = 0.25 rad, tdur = 10 sec, 
Ntot = 50, v' m — 1 keV, and p = 2.5 unless otherwise stated. We will see how the result 
changes if we change v' a , which ranges between 1 x 10~ 3 eV and 0.1 eV. We stress that 
the aim of this Letter is to show that the observed diversity of the prompt optical emission 
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can possibly be explained by our theoretical model. Therefore, a thorough study of the 
dependence of the other parameters is beyond the scope of this Letter and of future work. 

Before discussing further, one should note that since f(v') has a power-law form of 
f{v') oc u' 13 , Eqs. (CQ) and <^ provide us with 



F<? oc (T — t U) 



dep-> 



(6) 



which is the same form as a well known f ormula that is believed to describe the steep 



decay phase in the early X-ray afterglow (IKumar fc Panaitescu 



2000 



Lia ng et al. 



2006 



2006 



Butler fc Kocevski 



Nousek et al.ll2006l : lO'Brien et al.ll2006al : IZhang et al.ll2006l : lYamazaki et all 
20061 ). If the absorption frequency at the observer, v a [y a ~ 271/ if T ~ T s ( t a rt ), is below 
u Rl then (3 = —1/2 at ur, so that F^ j) rapidly decreases with r U) . On the other hand, if 

1 /2 

v a > v r-, then (3 = 5/2 at Ur, so that one can find oc r (j) . Keeping this in mind, we 
expect that the optical behavior for v a < vr is rather different from that for vr < u a , and 
that the observed diversity of prompt optical phenomena arises when v a varies across vr. In 
this Letter, we consider F Ux=w keV as a representative of the X-ray/7-ray light curve because 
as long as 1 keV < v < 100 keV, the light curve F v shows almost the same behavior. The 
following five cases (i)-(v), correspond to those in § [1] 

(i) Let us first consider the case 27^ < vr. In this case, we obtain the light curves in 
the R-band [vr = 1.7 eV) and the X-ray band [yx = 10 keV) as 



F (J) 



(7) 



'(j) 



(r s tart < T (j) < ^V'jvx) 

(7*4/ "x < t U) < r end ) 



where 



''"start 
''end 



(n 2 /r )(T, 



start tdep/ 

(c7 2 Ao)(T cn d-0 



= 1/2 , 
[1 + (7A0 t 



(9) 
(10) 



and we use the approximation 1 — /3 ~ l/2 / y 2 . Hence, both the optical and the X-ray light 
curves show rapid decay after the initial sudden rise, so that the optical behavior is similar 
to the X-ray behavior. When jV to t 3> 1, we see many spikes as observed for typical bursts. 
In particular, the time of the maximum flux is the same in the optical and the X-ray bands. 
In Figure la, we chose v' a = 1 x 10 -3 eV with other parameters being fiducial. The peak 
flux ratio is estimated as i? pea k = F^ k / F^^ ~ (^x/^r) 1 ^ 2 ~ 10 2 , which is larger than the 
observed value. If we consider the additional Band function component which dominates 
only in the observed X-ray and 7-ray bands, the value of -R pea k(~ several tens) can be closer 
to the observed one. 
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(ii) Next we consider the case vr < 2jv' a . Then we derive 

"* \ K^v'Jvr) 1 ' 2 ^ 2 (ju'J^r < r U) < r cnd ) ' 

while the X-ray light curve is again described by Eq. ([8]). The optical light curve has 
a gradually increasing part whose duration is ~ (l l 'a/ l/ R)( r o/ c 1 2 )- m particular, in the 
case of tdur ^ (l h 'a/ h 'R)(. r o/ c l 2 ), the difference between optical and X-ray behavior becomes 
significant. If iVtot ^ 1 and tdur < (7 Z7 a/ z/ fi)( r o/ c 7 2 ) ) all components are overlaid with 
each other at the time interval tdur < T < (•yb / ' a /b>R)(r Q /cy 2 ), resulting in a smooth light 
curve. On the other hand, since an X-ray light curve has short-duration, rapidly decaying 
pulses, we see less of a superposition effect than we do in the optical case. In Figure lb, we 
chose v' a = 0.1 eV with other fiducial parameters. One can see that the peak flux ratio is 
-Rpeak ~ 0(1), which is roughly consistent with the observation. 

(iii) In the intermediate case vr < 271/, we can see optical pulses synchronizing with 
the X-ray/7-ray pulses overlaid by a smooth component. In Figure lc, we chose v' a = 0.02 eV 
with other fiducial parameters. 

(iv) In the case of vr < 2^v' a , our model can also reproduce bright optical flash asso- 
ciated with e.g., GRB 990123. In Figure Id, we chose r = 3 x 10 15 cm, v' m — 1 eV, and 
u' a = 0.05 eV with other fiducial parameters. Since ji/ m <C v x = 10 keV, the spectrum at v x 
enters into the regime v' m < v' in Eq. ([3]), resulting in the smaller observed X-ray flux F vx 
compared with cases (i)-(iii). Then the large value of observed R pea k can be reproduced. 
Our result implies that exte rnal reverse shock emission is not necessarily required. Recently, 



Panaitescu fc Kumar! (120061 ) proposed a synchrotron-inverse Compton model to explain si- 
multaneously the optical flash and the 7-ray behavior of GRB 990123, which is similar claim 
as ours. Note that since u' m — 1 eV, an additional Band function component is necessary 
in order to match the observed spectrum in the range v > vx- Such a correction can be 
possible without changing the value of -R pea k ~ 10 2 . 

(v) Finally, if vr <C 27//, then the observed flux in the optical band becomes dim 
because of the steep slope below u a via the synchrotron self-absorption. 



3. Discussion 

In this Letter, assuming the synchrotron radiation with the absorption frequency v a 
around the R-band frequency ur, we have shown that the observed diversity of the prompt 
optical emission arises when v a varies and that our theoretical model explains the observa- 
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tion well. At present, the broad-band spectrum of the pr ompt emission of the GRB is un- 



Li & Sone 


2004; 


Zou et al. 


2006 



Panaitescu &: Kumarll2006l ). Our discussion can be generalized for such cases. If the comov- 
ing spectral function, f(v'), has the following form around a break frequency v' h ~ vr/^, 



(12) 



where s > 2 and q < 0, then, our conclusion derived in this Letter remains unchanged 
qualitatively. 

Our model will be tested by observing the spectral energy index, fto, in the optical band. 
When the optical light curve synchronizes with the X-ray/7-rays (i.e., ~ 27^ < ur), then 
the optical band is above the break frequency, u b , so that the energy index is flo ~ q < 0. 
On the other hand, if the optical light curve is rather smooth and does not synchronize 
with the X-ray/7-rays light curves, then the optical band is below v\> {vr < Vf,). Therefore, 
(3q ~ s > 2 is expected. Therefore, multicolor observation in the optical (and/or infrared) 
band is important to test our model. 

We would like to thank the anonymous referee and K. Ioka for useful comments. This 
work was supported by Grant-in-Aid for Scientific Research of the Japanese Ministry of 
Education, Culture, Sports, Science and Technology, 18740153 (R. Y.). 



REFERENCES 

Akerlof, C. et al. 1999, Nature, 398, 400 
Band, D. L. et al. 1993, ApJ, 413, 281 
Briggs, M. S. et al. 1999, ApJ, 524, 82 



Butler, N. R. & Kocevski, D. 2006, |astro-ph/ 0612564 

Fan, Y. Z. & Wei, D. M. 2004, ApJ, 615, L69 

Fan, Y. Z., Zhang, B., & Wei, D. M. 2005, ApJ, 628, L25 

Klotz, A. et al. 2006, A&A, 451, L39 

Kumar, P. & Panaitescu, A. 2000, ApJ, 541, L51 



- 7- 

Li, Z. & Song, L. M. 2004, ApJ, 608, L17 
Liang, E. W., et al. 2006, ApJ, 646, 351 

McMahon, E., Kumar, P., & Piran, T. 2006, MNRAS, 366, 575 

Meszaros, P. & Rees, M. 1999, MNRAS, 306, L39 

Nakar, E. & Piran, T. 2004, MNRAS, 353, 647 

Nousek, J. A. et al. 2006, ApJ, 642, 389 

O'Brien, P. T. et al. 2006a, ApJ, 647, 1213 

Panaitescu, A. & Meszaros, P. 2000, ApJ, 544, L17 

Panaitescu, A. & Kumar, P. 2006, |astro-ph/0612504 

Roming, P. W. A. et al. 2006, ApJ, 652, 1416 

Rykoff, E. S. et al. 2006, ApJ, 638, L5 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Sari, R. & Piran, T. 1999, ApJ, 520, 641 

Stanek, K. Z. et al. 2007, ApJ, 654, L21 

Toma, K., Yamazaki, R., & Nakamura, T. 2005a, ApJ, 620, 835 

Toma, K., Yamazaki, R., & Nakamura, T. 2005b, ApJ, 635, 481 

Vestrand, W. T. et al. 2005, Nature, 435, 178 

Vestrand, W. T. et al. 2006, Nature, 442, 172 

Wang, W. Y., Dai, Z. G., & Lu, T. 2000, MNRAS, 319, 1159 

Wei, D. M. 2007, MNRAS, 374, 525 

Yamazaki, R., Ioka, K., & Nakamura, T. 2004, ApJ, 607, L103 

Yamazaki, R., Toma, K., Ioka, K., & Nakamura, T. 2006, MNRAS, 369, 311 

Yost, S. A. et al. 2006, |astro-ph/0611414 

Zhang, B., Kobayashi, S., & Meszaros, P. 2003, ApJ, 595, 950 



Zhang, B. et al. 2006, ApJ, 642, 354 



Zhang, B. 2007, astro-ph/0701520 



Zou, Y. C, Xu, D., & Dai, Z. G. 2006, ApJ, 646, 1098 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 9- 




Fig. 1.— Examples of light curves. Solid and dotted lines are for the optical {vr = 
1.7 eV) and X/7-ray {u x = 10 keV) bands, respectively. The adopted parameters are 
(a) i/ a = lx 1(T 3 eV, (b) 1/ = 0.1 eV, (c) v' a = 0.02 eV, and (d) v' a = 0.05 eV, v' m = 1 eV, 
r = 3 x 10 15 cm, with other fiducial parameters (r = 1 x 10 15 cm, 7 = 100, A6> to t = 0.25 rad, 
tdur = 10 sec, iVtot = 50, v' m — \ keV, and p = 2.5). 



